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Abstract 

The effect of plasma screening on the electronic and vibrational properties of the Ha^ molecular 
ion was analyzed within the Born-Oppenheimer approximation. When a molecule is embedded 
in a plasma, the plasma screens the electrostatic interactions. This screening is accounted for in 
the Schrodinger equation by replacing the Coulomb potentials with Yukawa potentials that 
incorporate the Debye length as a screening parameter. Variational expansions in confocal 
elliptical coordinates were used to calculate energies of the l^cr^. and the ^pa^ states over a 

range of Debye lengths and bond distances. When the Debye length is comparable to the 
equilibrium bond distance, the plasma screening reshapes the potential energy curve. Expectation 
values, dipole polarizabilities and spectroscopic constants were calculated for the Isa state. 
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1. Introduction 

The interaction of atomic systems with plasmas has been studied extensively [1-9] and that 
work has recently been extended to molecular systems [10, 11]. Plasma embedding is modeled by 
replacing the Coulomb potentials of the free system with Yukawa potentials [1,2] where the 
screening of the plasma is parameterized by the Deb ye length. It has been shown [10-14] that 
decreasing the Debye length is accompanied by a decrease in the binding energy. In this work we 
investigated the effect of plasma embedding on the potential energy curves of the Iscr^ and 2o(X, 

states of the H2* molecular ion. was chosen because atomic and molecular hydrogen are 
abundant in interstellar matter and are of interest in astrochemistry[15]. The Isa state is the 

prototype for molecular bonding while the ^pa^ state is an example of a van der Waals molecule 

[16] that is bound at internuclear separations greater thani? = 10.41 ao- The effects of plasma 
embedding on the bonding behavior of these two states are described in this Letter. We also 
investigated the effects of Debye screening on the quadrupole moment integrals, parallel and 
perpendicular dipole polarizabilities, Dunham coefficients, harmonic force constants and 
harmonic fi-equencies of the Is a state. 



2. Calculational method 

Within the Bom-Oppenheimer approximation, the Schrodinger equation in atomic units for 
H2^is 
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where and are the distances from the electron to nuclei A and B and R is the internuclear 
separation. This equation is most easily solved in confocal elliptical coordinates defined as 
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where ^is the angle of rotation about the internuclear axis. The Schrodinger equation for cr states 
of is now 
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Equation 3 is separable and has an analj'tic solution as the product of two infinite series [17,18]. 

When an atom or molecule is embedded in a hot, dense plasma, the electrostatic interactions 
between the particles are screened by the plasma. The Debye-Hiickel theory - first formulated in 

the theory of electrolytes [19] - is widely used for modeling plasma screening because the 
required integrals are analytic. This Debye screening is included in the Schrodinger equation by 
replacing the Coulomb potentials by Yukawa potentials [1], 

r 

V{r)^—, (4) 
r 

where D, the Debye length, is a measure of the distance over which the Coulomb potential is 



killed off by the polarization of the plasma. Physically/) oc -J¥7n , where T is the temperature of 
the plasma and n is its number density. In the limit of large D the Yukawa potential approaches 
the Coulomb potential and the plasma-embedded system approaches fi-ee H2^. However, if D is 
comparable to the dimensions of the system, Debye screening significantly affects the 
interactions of the particles. 

Incorporation of the Yukawa potential into the Schrodinger equation requires replacement of 
the potential energy portion of equation (3) by 



V = -- 
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The resulting Schrodinger equation is non-separable and must be solved approximately. 
The present work used a variational treatment with a trial wavefunction 

N m 



(6) 



where the P^^^^ are the Legendre polynomials of degree I. This form of wavefunction was 

chosen so as to include terms similar to those used in the exact calculation [17,18]. For states of 
gerade parity, I was restricted to even values while for states of ungerade parity, i was odd. 
The Cf were determined by solving the secular equation and or was hand optimized to give ten 

digit convergence at each value of/? and D. Using these wavefunctions to calculate the ground 
and excited state energies for unscreened gave variational energies that agreed with the exact 
energies [18] through 1x10"^ £'h. 
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For each value of Z) we fitted a polynomial through the points around the minimum in the 
energy curves to find the equilibrium separation, R^, and the energy at R^. 
The quadrupole moment integrals 

(z^) = (v^|z», (x^) = (5^|x», (7) 

depend on the wavefunction at distances from the nuclei that are greater than those which 
effectively determine the energy [20] and their dependence on D is thus of some interest. They 
can also be used in the Kirkwood approximation [21] 

4{z'f<a^^,4(x'f<a^ (8) 

to obtain lower bounds to a^^ and , the parallel and perpendicular dipole polarizabilites. 
The polarizabilities can be obtained from perturbation theory as 

a^^=-2{i^\z\i^^^),a^=-2{i^\x\i^-^), (9) 

where and y/^ are the first-order corrections to y/ under the perturbation operators z and x, 
respectively. Hylleraas variational perturbation theory [22] has been shown [23] to be an 
accurate method to find i^^^ and y/^ and lower bounds to ay and for H2^. 

The Dunham parameterization [24] of the potential energy curve corresponds to a fourth- 
degree pol5Tiomial fit around the minimum in the potential energy where 

E{R) = E{R^) + A,p'(l + A,p + A,p') (10) 

and 

p = ^. (11) 

Re 

The fitting parameters can then be used to find the harmonic force constant (^e) and the harmonic 
frequency (coe). 

3. Results and discussion 

ha^ energies were calculated over the range i? = 1 - 10 ao for-D = 50, 20, 10, 5, 2 and 1. 

The calculation used the metric i + m < 10, m = even, giving a 36-term waveftmction. Iscr^ 

energies are given in Table 1 and are shown graphically in Figure 1. The effect of screening is 
more easily seen by calculating the potential energy curve U {D,R^ where 
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U{D,R) = E^^ {D,R)-E^ (D) , (12) 

and gives the energy of the screened molecular ion relative to (£)) , the energy of a screened 
hydrogen atom with the same value of Z) . The depth of these curves at Re is the dissociation 
energy of the embedded molecular ion, D^. E^ (Z)) was calculated using the p-FEM method 

[25]. Results from the p-FEM method were consistent with [14] and with the recent work of Paul 
and Ho [26]. 

Potential energy curves for Z> = oo, 5, 2 and 1 are shown in Figure 2. The curves for D = 50 
and D = 20 were not shown as they are so close together at this scale as to be indistinguishable 
from the D = co curve. For D>5ao, the effect of screening is small, with the change in binding 
energy of the molecular ion very nearly equal to the change in the energy of the Yukawa- 
screened hydrogen. For D = 2 and 1, the potential energy curves become markedly shallower and 
the minima in the curves shift to larger R. The changes in the potential energy curves result from 
an increase in the elecfronic energies (which include the kinetic energy, the elecfron-nuclear 
attraction and the energy of the screened hydrogen atom) offset by a slightly smaller decrease in 
the intemuclear repulsion term, resulting in an overall decrease in stability. This decrease in 
stability with decreasing D is consistent with the results observed by Mukherjee et al. [10]. 

2pa^ energies were calculated over the range /? = 10 - 20 ao for Z) = 50, 20, 10, 5, 4, 3 and 

2. For Z) = 2, the excited state was non-bonding relative to the Yukawa-screened hydrogen atom. 
The metric was increased to 12 with m = odd in order to get convergence at small D. This gave a 
49-term wavefunction. 

2 pa^ energies are given in Table 2 and the potential energy curves for D = co, 10, 5, 4, 3 

and 2 are shown in Figure 3. The 2 pa^ state is very loosely bound even at large D and bonding 

occurs only at intemuclear separations greater than 10 uq. For Z> > 10. the dissociation energy 
initially increases and the equilibrium intemuclear separation decreases as the Debye screening 
increases. For D< 10, the dissociation energy decreases accompanied by an increase in the 
intemuclear separation. 

As D decreases in the region where Z) > 10, the intemuclear repulsion term decreases more 
rapidly than the elecfronic term increases, resulting a small increase in stability. For D<10, the 
2 pcr^ state behaves like the Iscr^ state. In both cases, the changes in the dissociation energy are 
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of the order of 10""* £'h while the equilibrium intemuclear separation varies from 12.08 to 13.90 

Values of the equilibrium intemuclear separation and dissociation energy are included in 

Tables 1 and 2. 

Tables 3-5 show the dependence on Z) of a variety of Isa^ electronic and spectroscopic 

properties. In all cases, they were evaluated at the value of corresponding to the equilibrium 
nuclear separation for that choice of D. Table 3 shows selected expectation values as a function 

ofD. For > 5 ao, {V) , (z') and ( x^^ are relatively constant. However, when D becomes 
comparable to the intemuclear separation (7? « 2 ao), (V) drops significantly, while ^z^^ and 

i^x^^ show a marked increase. Since ^z^^ and ^x^^ measure the extent of the electron 

distribution and, through Eq. (8), provide lower bounds to the polarizability, we interpret this 
increase as consistent with a more diffuse, and thus more polarizable, electron distribution. 

Table 4 shows the parallel and perpendicular dipole polarizabilities along with the Kirkwood 
lower bounds of Eq. (8). We have also shown the average polarizability 



and the anisotropy 



« = ^(«ll+2«x), (13) 



a , CKj^ and a show the same trends as ^z^ ^ and ^x^ ^ , consistent with the conclusion of Kar and 

Ho [1 1] that the system becomes more polarizable as D decreases. The anisotropy decreases with 
with decreasing D and the electron distribution becomes more diffuse, approaching a near- 
spherical geometry for small D. 

The Dunham parameters Aq, Ai and A2 are shown in Table 5 along with the harmonic force 

constant (^e = 2AqIR^) and the harmonic frequency {^(0^= ^^/3.633xlO~^* j . Again we see a 
significant decrease in and cot at small D. Our free system values for Aq = 0.2054 Eh, A\ = - 
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1.7571 and ^2=2.1 264 for Z> = oo are in good agreement with the values of [27] (Aq = 0.2053 £"11, 
^1 = -1.7363 and ^2 =2.1329). 

In summary, the electronic and vibrational properties of embedded in a Debye plasma 
are significantly affected when D » R^. The dissociation energy is reduced, the equilibrium 
intemuclear separation increases, the polarizability increases and vibrational force constant is 
reduced. While the dipole polarizabilites calculated in this work were for a static electric field, it 
would be of interest to investigate effect of plasma embedding on the dynamic polarizabilies. 
This investigation is one of our goals. 
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^ics 111 Hi lui lL^ \j g 


embedded in a Debye plasma 












D — CO 


Z) = 50 


£) = 20 


r\ 1 A 

U — \\} 


L> — J 


U — 1 


n — 1 


1 


-0.45178631 


-0.43203847 


-0.40333621 


-0.35782414 


-0.27479473 


-0.07927245 


0.10997427 


2 


-0.60263421 


-0.58289414 


-0.55422709 


-0.50881491 


-0.42606393 


-0.23172997 


-0.04520628 


3 


-0.57756286 


-0.55783451 


-0.52922719 


-0.48402329 


-0.40210546 


-0.21366909 


-0.04422103 


4 


-0.54608488 


-U. JZOjO / / / 


-U.HV / i 07 J 


-0.45282386 


-0.37175069 


-0.18892452 


-0.03439361 


5 


-0.52442030 


-0.50471245 


-0.47621058 


-0.43137647 


-0.35090744 


-0.17165594 


-0.02705181 


6 


-0.51196905 


-0.49226768 


-0.46379661 


-0.41905766 


-0.33889097 


-0.16116181 


-0.02221179 


7 


-0.50559400 


-0.48589640 


-0.45744133 


-0.41274289 


-0.33266080 


-0.15518356 


-0.01901731 


g 


-0 50257039 


-0.48287454 


-0.45442560 


-0 40973585 


-0 32963609 


-0 1 51 90162 


-0 01684549 


9 


-0.50119545 


-0.48150009 


-0.45305208 


-0.40835735 


-0.32821205 


-0.15013570 


-0.01531628 


10 


-0.50057873 


-0.48088349 


-0.45243369 


-0.40772991 


-0.32754236 


-0.14919346 


-0.01420328 


00 


-0.50000000 


-0.48029611 


-0.45181643 


-0.40705803 


-0.32680851 


-0.14811702 


-0.01028579 




1.996825 


1.996933 


1.997465 


1.999247 


2.006098 


2.060332 


2.360900 




-0.60263461 


-0.58289450 


-0.55422730 


-0.50881482 


-0.42606567 


-0.23188935 


-0.04831233 


D, 


0.10263461 


0.10259839 


0.10241087 


0.10175678 


0.09925716 


0.08377233 


0.03802654 




Fig 1. Is a energies 
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Fig 2, ha^ U{D,R) for Z) = oo, 5, 2 and 1 . 



Table 2. Energies in for 2p (t, 


H^' embedded in a Debye plasma 












R («o) 


D = X 


Z5 =50 


D =20 


D = 10 


£1=5 


£1=4 


£1=3 


£1=2 


10 


-0.49990107 


-0.48020631 


-0.45175695 


-0.40704313 


-0.32677122 


-0.29081730 


-0.23656384 


-0.14732128 


11 


-0.50002442 


-0.48032909 


-0.45187710 


-0.40715783 


-0.32689160 


-0.29095505 


-0.23675433 


-0.14770032 


12 


-0.50005789 


-0.48036205 


-0.45190768 


-0.40718307 


-0.32691837 


-0.29099315 


-0.23682727 


-0.14789715 


13 


-0.50005947 


-0.48036316 


-0.45190655 


-0.40717685 


-0.32691160 


-0.29099390 


-0.23685097 


-0.14800024 


14 


-0.50005150 


-0.48035475 


-0.45189606 


-0.40716168 


-0.32689521 


-0.29098278 


-0.23685521 


-0.14805475 


15 


-0.50004206 


-0.48034488 


-0.45188426 


-0.40714571 


-0.32687813 


-0.29096965 


-0.23685261 


-0.14808364 


16 


-0.50003368 


-0.48033606 


-0.45187369 


-0.40713147 


-0.32686322 


-0.29095788 


-0.23684826 


-0.14809905 


17 


-0.50002688 


-0.48032881 


-0.45186483 


-0.40711946 


-0.32685102 


-0.29094825 


-0.23684401 


-0.14810730 


18 


-0.50002155 


-0.48032295 


-0.45185750 


-0.40710944 


-0.32684124 


-0.29094063 


-0.23684034 


-0.14811172 


19 


-0.50001742 


-0.48031813 


-0.45185135 


-0.40710099 


-0.32683336 


-0.29093456 


-0.23683717 


-0.14811407 


20 


-0.50001420 


-0.48031396 


-0.45184597 


-0.40709367 


-0.32682679 


-0.29092949 


-0.23683423 


-0.14811528 


00 


-0.50000000 


-0.48029611 


-0.45181643 


-0.40705803 


-0.32680851 


-0.29091959 


-0.23683267 


-0.14811702 


(«0) 


12.422056 


12.396172 


12.284674 


12.081375 


12.082440 


12.392866 


13.901264 




E(R,) 


-0.500061 


-0.480365 


-0.451910 


-0.407184 


-0.326919 


-0.290996 


-0.236855 






-0.000061 


-0.000069 


-0.000093 


-0.000126 


-0.0001 1 1 


-0.000077 


-0.000023 
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Fig 3. 2 pcr„ [/(Z),i?) for Z) = oo, 10, 5, 4, 3 and 2. 



Table 3. Isag 


expectation values for potential 


energy and quadrupole integrals 


D 


<V> 


<z^> 


<x'> 


00 


-1.205268 


1.108927 


0.641143 


50 


-1.185273 


1.109463 


0.641506 


20 


-1.155327 


1.112138 


0.643330 


10 


-1.105631 


1.121110 


0.649482 


5 


-1.007629 


1.153606 


0.672235 


2 


-0.728481 


1.371570 


0.831278 


1 


-0.318519 


2.804797 


1.996692 
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Table 4. 1* ct^ dipole polarizabilities and Kirkwood lower boimds (a o^) 



D 




«-Li.b. 




a 


A* 


00 


4.918876 5.057482 


1.644256 


1.754475 


2.855477 


0.385576 


50 


4.923633 5.062401 


1.646120 


1.756608 


2.858539 


0.385487 


20 


4.947400 5.086966 


1.655496 


1.767330 


2.873875 


0.385036 


10 


5.027550 5.169825 


1.687306 


1.803680 


2.925728 


0.383511 


5 


5.323223 5.475791 


1.807601 


1.941421 


3.119544 


0.377659 


2 


7.524813 7.776257 


2.764093 


3.060827 


4.632637 


0.339291 


1 


31.467546 35.103927 


15.947112 


20.734794 


25.524505 


0.187651 


Tabic 5. l.va^, Dunham parameters, 


harmonic force 


constant (/v\,) and harmonic frequency (co^) 


D 








^e(^h/ao^) 


CO ^(cm^) 


00 


0.2054 


-1.7571 


2.1264 


0.1030 


2325 


50 


0.2054 


1.7573 


2.1273 


0.1030 


2324 


20 


0.2052 


-1.7580 


2.1316 


0.1028 


2322 


10 


0.2044 


-1.7608 


2.1467 


0.1022 


2316 


5 


0.2012 


-1.7725 


2.2085 


0.1000 


2290 


2 


0.1761 


-1.8435 


2.7753 


0.0831 


2088 


1 


0.0936 


-2.1862 


3.7030 


0.0338 


1332 
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